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Abstract 

The measurements of the fluxes of atmospheric neutrinos give evidence for the 
disappearance of muon neutrinos. The determination of the dependence of the 
disappearance probability on the neutrino energy and trajectory allows in princi- 
ple to establish unambiguosly the existence of neutrino oscillations. Alternative 
mechanisms for the disappearance of the neutrinos have been proposed, but do not 
provide a viable description of the data, if one includes both events where the neu- 
trinos interact in the detector and ^-induced upward going muons. The proposed 
mechanisms differ in the energy dependence of the disappearance probability and 
the upward going muon data that are produced by high energy neutrinos give a 
i-^ ■ crucial constraint. 

> : 

•rH , 

^ : 1 Introduction 

The measurements of the fluxes of atmospheric neutrinos by the Super-Kamiokande (SK) 
experiment |], |2|, || [|, || show evidence for the disappearance of muon (anti)-neutrinos. 
The same indication comes from the older data of the Kamiokande || and 1MB |7] ex- 
periments and the recent ones of Soudan-2 || . Also the results recently presented by the 
MACRO collaboration || [TD| indicate a suppression of the muon (anti)-neutrino flux. 

The simplest explanation of the data is the existence of <-> v T oscillations . In the 
framework of flavor oscillations one should consider the more general case of three flavors 



TTf| (with the Chooz experiment [12| giving important contraints to the electron neutrino 



transitions), and could also envisage more complex scenarios involving sterile states |T3 
We will not pursue these possibilities here, and we will adopt instead the simplest scenario 
of two-flavors oscillations as a prototype model that, as we will see, is able to describe 
successfully the experimental data. 

We will instead investigate if other forms of 'new physics' beyond the standard model, 
different from standard flavor oscillations, can also provide a satisfactory description of 
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the existing data. Indeed several other physical mechanisms have been proposed in the 
literature as viable explanations of the atmospheric neutrino data. In this work we will 
consider three of these models: neutrino decay [14 
(FCNC) || 



16fl , and violations of the equivalence principle [[3j], [18 



flavor changing neutral currents 
equivalent ly, 



or, 



of Lorentz invariance jl9" |. All these model have the common feature of 'disappearing' 



muon neutrinos, however the probability depends in different ways on the neutrino energy 
and path. To discriminate between these models a detailed study of the disappearance 
probability P and of its functional form is needed. 

In this work in contrast with previous analyses we will argue that the present data 
allow to exclude the three 'exotic' models, at least in their simplest form, as explanations 
of the atmospheric neutrino problem. This is mainly due to the difficulty that these 
models have to fit at the same time the SK data for leptons generated inside the detector 
(sub- and multi-GeV) and for up-going muons generated in the rock below it. 



2 Data 

In fig. 1, 2 and 3 we show (as data points with statistical error bars) the ratios between 
the SK data and their Montecarlo predictions calculated in the absence of oscillations 
or other form of 'new physics' beyond the standard model. In fig. 1 we show the data 
for the e-like contained events, in fig. 2 for /i-like events produced in the detector, and 
in fig. 3 for upward-going muon events, as a function of zenith angle of the detected 
lepton. In each figure we include four lines: the dotted line (a constant of level unity) 
corresponds to exact agreement between data and no-oscillation Montecarlo, including 
the absolute normalization. The dot-dashed lines correspond to the assumption that there 
is no deformation in the shape of the zenith angle distributions, but that one is allowed 
to change the normalization of each data sample independently. The values obtained are: 
1.16 for e-like sub-GeV, 1.21 for e-like multi-GeV, 0.72 for /i-like sub-GeV, 0.74 for /i- 
like multi-GeV, 0.56 for stopping upward-going muons, and 0.92 for passing upward-going 
muons. For two sets of data (sub-GeV and multi-GeV /i-like events) the constant shape 
fits give very poor descriptions (x 2 = 26 for the sub-GeV and 33 for the multi-GeV for 4 
d.o.f). Also the zenith angle shape of the passing upward-going muons is not well fitted 
by the no-oscillation Montecarlo (x 2 = 17 for 9 degrees of freedom). The electron data do 
not show clear evidence of deformations, although the constant shape fit for the sub-GeV 
events (x 2 = 9.7 for 4 d.o.f.) is rather poor. 

The normalizations of the different data sets are of course strongly correlated, and 
therefore it is not reasonable to let them vary independently. The other extreme option, 
that we will adopt in this work for simplicity, is to use one and the same parameter to fix 
the normalization of the six data samples. The result for constant shapes (i.e. assuming 
no 'new physics' beyond the standard model) is represented by the dashed lines in fig. 1, 
2 and 3 corresponding to a value 0.84 and a very poor x 2 — 280 for 34 d.o.f.). 

The full lines in the figures correspond to our best fit assuming z/ M <-> v T oscillations 



2 



with maximal mixing. We define the \ 2 as follows: 



X 



E 



Nj - a iVf (N?$ c /N* 



-i 2 



j,0J 



(1) 



In (P the summation runs over all data bins, Nj is the SK result for the j—th. bin, cxj 
its statistical error, Nj h our prediction, N* h our prediction in the absence of oscillations, 
N?mc the no-oscillation prediction of Super-Kamiokande, and a allows for variations 
in the absolute normalization of the prediction. We have rescaled our prediction to the 
SK Montecarlo because we do not have a sufficiently detailed knowledge of the detector 
response (e.g. number of detected rings) and efficiency. For the same same input neutrino 
spectra the difference between our no-oscillation calculation (see pO for a description) 
and the SK Montecarlo result is approximately 10%, 

For our best fit the values of the relevant parameters are a = 1.15 and Am 2 = 
3.2 1(T 3 eV 2 . The X 2 is 33.3 for 33 d.o.f. 

Our definition of the x 2 is somewhat simplistic. We do not take into account the 
contribution of systematic errors, either in the data or in the theory. The assumption 
of a common a for e-like and //-like events corresponding to different energy regions is 
certainly too strict. It is therefore remarkable that this fit is so good, and essentially in 
agreement (same normalization and very near Am 2 value) with the much more elaborate 
fit in 0. 

In the rest of this paper we will consider other, 'exotic' models and we will find that 
they are not able to provide a satisfactory fit to the same data. 



3 Models 



We briefly recall the essential points of the models we are discussing. 

For the usual, two-neutrino flavor oscillations the 'disappearance probability' P is 
given by: 

"Am 2 L 



pc 



sin 29 sin 



E v 



(2) 



with the very characteristic sinusoidal dependence on the ratio L/E v . 

In the simplest realization of neutrino decay, neglecting the possibility of the simulta- 
neous existence of neutrino oscillations, the disappearance probability is given by: 



P 



dec 



exp 



m. 



L 



T„ E,. 



(3) 



still depending on the ratio between neutrino pathlength and energy L/E u , but with a 
functional form different from (@). 

If flavor changing neutral currents contribute to the interaction of neutrinos with or- 
dinary matter, a non trivial flavor evolution will develop even for massless neutrinos as 



3 



originally noted by Wolfenstein |]2T|. There are several theoretical models generically pre- 
dicting nondiagonal neutrino interactions with matter. In particular such models have 
been proposed as a possible consequence of i?-parity violating interactions in supersym- 
metric models and suggested as solutions of both the solar |22|, 15] and atmospheric | 15| , |16 
neutrino problems. Let us call V a p the effective potential that arises from the forward scat- 
tering amplitude of a neutrino with a fermion /: v a + / — > i/p + /. In the standard model 



V, 



V Tfl = 0, and = V TT = V2G F T 3 (f L ) 



Nf where Gp is the Fermi constant, Nf 



is 



the number density of the fermion / and T^fa) is the third component of the fermion's 
weak isospin. Since the effective potentials for muon and tau neutrinos are identical, there 
is no effect on standard oscillations. However, if the scattering amplitudes are different 
from those predicted by the standard model, and if flavor changing scattering can occur, 
then the effective potential acquires non diagonal terms V^ T = V Ttl = v^G^eiV/, and 
different diagonal elements (with V TT — = y^GVe'iV/), and there will be a nontrivial 
flavor transition probability even for massless neutrinos. After the crossing of a layer of 
matter with a column density 



dL> Nf(L') 



the transition probability is: 



P = P, 



FCNC 



4e 2 



4e 2 + 
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(4) 



(5) 



The probability has again an oscillatory form, however in this case the role of L/E u is 
taken by the column density Xf and there is no dependence on the neutrino energy. 

If the gravitational coupling of neutrinos are flavor dependent (implying a violation of 
the equivalence principle) mixing will take place for neutrinos traveling in a gravitational 
field even for massless neutrinos 
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The neutrino states with well defined coupling to 
the gravitational field define a 'gravitational basis' related to the flavor basis by a unitary 
transformation. The effective interaction energy matrix of neutrinos in a gravitational 
field can be written in an arbitrary basis as 



H = -2\<P(r)\ E v (l + f) 



(6) 



where E v is the neutrino energy, cf){r) = — \<fi(r)\ is the gravitational potential, and / is 
a (small, traceless) matrix that parametrize the possibility of non-standard coupling of 
neutrinos to gravity and is diagonal in the gravitational basis. 

Much in the same way as in the previous cases, the noncoincidence of gravitational 
and flavor eigenstates determines mixing and flavor transitions. Considering the simple 
case of two flavors and assuming a constant gravitational potential \<p\, the transition 
probability takes the form 



p _ pgrav 



sin 2 (26 



g) sin 



E V L\. 



(7) 



where Qq is the mixing angle and 5 is the difference between the coupling to gravity of the 
gravitational eigenstates. Note that in this case the argument of the oscillatory function 
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is proportional to the product of the neutrino energy and pathlength, whereas for the 
standard flavor oscillations it is the ratio of the same quantities that matters. 

Equations (0), (|3]), (|5|) and (0) are the disappearance probabilities for the four mech- 
anisms that we will confront with the experimental data. 

4 Flavor Oscillations 

It is interesting to discuss how the usual flavor oscillations can successfully reproduce the 
pattern of suppression measured for the different event samples. The events detected in 
one particular bin are produced by neutrinos with a predictable distribution of E v and 
pathlength L, and therefore of L/E v , the significant quantity in flavor oscillations. In 
fig. 4, in the top panel we show as a function oiL/E v the survival probability correspond- 
ing to maximal mixing and Am 2 = 3.2 x 10~ 3 eV 2 (our best fit point). Also shown with a 
dashed line is the survival probability for neutrino decay that we will discuss in the next 
section. 

In the second panel we show the L/E u distributions of sub-GeV /x-like events in the 
five zenith angle bins used by the SK collaboration: cos 9^ G [—1,-0.6], [—0.6,-0.2] 
[—0.2,-1-0.2], [0.2,0.6] and [0.6,1.0] (corresponding to the thick solid, thick dashed, thin 
dot-dashed, thin dashed and thin solid line). 

In the third panel we show the corresponding distributions for multi-GeV /i-like events 
(same coding for the lines). 

In the fourth panel we show the L/E v distributions for upward going muons that stop 
in the detector in the zenith angle bins: cos# M G [—1,-0.8], [—0.8,-0.6] [—0.6,-0.4], 
[—0.4, —0.2] and [—0.2, 0.0] with the corresponding lines ordered from right (higher values 
of L/E v ) to left (lower values of L/E v ). 

In the last panel we show the same distributions for passing upward going muons in 
ten zenith angle bins, cos# M G [—1, —0.9], . . ., [—0.1,0.0]. 

Some remarks can be useful for an understanding of the distributions shown in fig. 4. 
For the sub-GeV events, one can see that the parent neutrinos have L/E u spread over a 
broad range of values. This is due to the poor correlation between the neutrino and muon 
directions (9 VI1 ) ~ 53°. 

For multi-GeV data the distributions are much narrower, reflecting the tighter corre- 
lation between the neutrino and muon directions, {9 V ^) ~ 13°. Note also that the peaks 
in the LjE v distributions corresponding to sub-GeV and multi-GeV events in the same 
zenith angle interval are at slightly different points because of the different energy of the 
parent neutrinos. 

For up-going stopping muons the width of the distribution is wider than in the multi- 
GeV case. The correlation between the muon and neutrino directions (9 VfJ ) ~ 10° is 
actually better, but the width of the distribution reflects the wider energy range of the 
neutrinos contributing to this signal. Passing muons are nearly collinear with the parent 
neutrinos ((9 Ufl ) ~ 2.9°), but the large energy range of the neutrinos that extend over 
nearly two decades (E v ~ 10-10 3 GeV) results in a wide L/E v distribution. 
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All curves in the lower four panels of fig. 4 are normalized to unit area. In order 
to obtain the suppression due to oscillations in a particular bin, one has to perform the 
integral: 

r dN° 

iVf c (sm 2 29, Am 2 ) = J dx [1 - P osc (x, sin 2 29, Am 2 )] (8) 

Comparing the survival probability with the L/E v distributions it is easy to gain a 
qualitative understanding of the effects produced. For Am 2 ~ 3 • 10~ 3 eV 2 , neutrinos 
with L/E u < 10 2 Km/GeV have a survival probability close to unity and do not oscillate, 
while for neutrinos with L/E v > 10 3 Km/GeV, averaging over the rapid oscillations, 
the survival probability becomes one half for maximal mixing. We recall that horizontal 
neutrinos travel an average pathlength of ~ 600 Km. 

Taking into account the L/E u distributions of the different set of events one can see 
that all zenith angle bins of the muon sub-GeV events are somewhat suppressed, because 
even vertically downward going muons can be produced by upgoing neutrinos. 

For multi-GeV events, with the tighter correlation between the neutrino and muon 
directions, the two up-going bins are suppressed by the 'average' factor ~ 0.5, the two 
down-going bins are left unchanged and the horizontal muons have an intermediate sup- 
pression. 

The up-going stopping muons are always suppressed by a factor ~ 1/2, except for the 
bin nearest to the horizontal. 

For the up-going passing muons the larger average energy and therefore smaller L/E u 
explains the smaller suppression and its pattern, varying from nearly unity for the hori- 
zontal bin to a maximum of ~ 0.65 for the vertical one. 



5 Exotic Models 
5.1 Neutrino Decay 

Fitting the sub-GeV and multi-GeV data of Super-Kamiokande with the simplified model 
of muon neutrino decay (that neglects mixing) given in ([3]), we find a minimum in the 
X 2 for a value T v jm v = 8900 Km/GeV (with a = 1.07). This is in good agreement 
with the results of ([TJj]. The authors of this reference have as a best fit point T v /m v ~ 
12800 Km/GeV, with a small mixing angle sin 2 29 ~ 0.06. The curve describing the decay 
probability for our best fit is shown as the dashed line in the top panel of fig. 4. 

It is simple to have a qualitative understanding of the value of r v /m v that provides 
the best fit. One needs to suppress by a factor ~ 0.5 the up-going multi-GeV muons that 
have (L/E v ) ~ 10 35 Km/GeV (see fig. 4). 

The inclusion of decay results in x 2 — 71 (for 18 d.o.f.), a very significant improvement 
over the value 234 (19 d.o.f.) of the "standard model", but still significantly worse than 
the value x 2 — 25 of the <-> v T flavor oscillation fit to the same set of data. 

For a value of r v /m v of the order of what is given by our fit to the sub-GeV and 
multi-GeV data, one expects a much smaller suppression of the high energy passing up- 
going muons (as already noted in [14]). In fact including also the 15 data points of the 



6 



up-going muons in a new fit, the best fit point becomes r v jm v = 10000 Km/GeV (similar 
to the previous one), but x 2 increases to the much higher value 140. 

5.2 Violation of the equivalence principle 

Performing a fit to the sub-GeV and multi-GeV data of Super Kamiokande with the 
disappearance probability given by (|7]) and with maximal mixing (9q = vr/4), we find a 
minimum in the \ 2 for a value $\<$>\ = 4 • 10~ 4 Km _1 GeV _1 (with a = 1.10). The x 2 f° r 
this fit is 35 for 18 d.o.f., still a very significant improvement over the standard model 
case, but not as good as the flavor oscillations result. The survival probability given by 
our best fit is shown in the top panel of fig. 5. 

The reason of the poor x 2 can qualitatively be understood looking at fig. 5. This figure 
is the equivalent of fig. 4, in the sense that the four lower panels show the distributions in 
the variable that is relevant in this case, namely L ■ E u . The distributions in this variable 
for the sub-GeV and multi-GeV events have shapes similar to the corresponding ones 
in L/E u , because the width of the distributions is mostly determined by the spread in 
pathlength L. However the average value of the L ■ E v of the sub-GeV events is lower 
than the corresponding one (same zenith angle bin) for multi-GeV events, the opposite 
of what happens in the L/E u distributions, see fig. 4. Therefore, parameters describing 
well multi-GeV events will generally produce too low a suppression for sub-GeV events 
or viceversa. 

It can be argued (as the authors of reference |24| do) that taking into account sys- 
tematic uncertainties the model defined by equation (^) provides a good fit to the data, 
however this is not the case if upward-going muons are included in the picture. This 
should be evident looking at the lower panels in fig. 5. Upward-going muons are pro- 
duced by high energy neutrinos and the frequent oscillations do imply a suppression by 
50% of passing (and stopping) muons, with no deformation of the zenith angle distribu- 
tion. This is in disagreement with the corresponding data: in fact, trying to fit all the 
data together we obtain similar best fit parameters, 8\<f>\ = 4.5 • 10 -4 Km _1 GeV _1 and 
a = 1.145, but with a very bad x 2 = 142.7 for 32 d.o.f. (the contribution of passing 
upward-going muon data being ~ 100). 

5.3 Flavor Changing Neutral Currents 

In the case of neutrino transitions produced by flavor changing neutral currents, the role 
oiL/Ey is replaced by X, the column density. This has the fundamental consequence that 
there is no energy dependence of the flavor conversion. Moreover since air has a density 
much lower that the Earth's, the transitions do not develop during the neutrino path in 
the atmosphere, and therefore down-going neutrinos are unaffected. Note also that there 
is not a simple relation between the zenith angle 9 V and the pathlength L because of 
fluctuations in the neutrino birth position. However, due the air low density, the zenith 
angle 9 U does define the column density X with a negligible error: the entire down-going 
hemisphere corresponds to X ~ and to a vanishing transition probability. 
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Performing, as before, a fit to the sub-GeV and multi-GeV data of Super Kamiokande 
with the disappearance probability given by (j|) and assuming scattering off down quarks 
and e' — (that is maximal mixing), we obtain a best fit value e = 0.4 and a = 1.08 
corresponding to a minimum x 2 = 38. With increasing e the oscillations become more 
frequent, and essentially all values e > 0.4 give comparable fits, since for these large values 
the oscillations can be considered as averaged in the entire up-going hemisphere. 

The authors of reference fll~6| , exploring the parameter space (e, e') find two solutions: 
(a): (0.98,0.02) and (b): (0.08,0.07), that are plotted in the upper panels of fig. 6. The 
first solution corresponds to the one that we have found, considering the slow variation 
of x 2 with e in the large e region. The x 2 found by the authors of Jl6[ is however better 
that what we find, indeed as good as in the flavor oscillation model. 

We do find that fitting the muon data only, without considering the constraint on 
the normalization coming from the electron data, the FCNC model gives an excellent 
fit, indeed as good or better than the flavor oscillation model. The reason why, in our 
fitting procedure, the FCNC model gives not as good a fit originates from the fact that 
the theoretical average value of the suppression for both sub-GeV and multi-GeV muon 
events for the best fit parameters is ~ 0.75, corresponding to no suppression in the down- 
going hemisphere and ~ 0.5 in the opposite one. The data || [3| for the double ratio 
R = (ji/e) Data /(ji/e) M c- Rsub = 0.61 ± 0.03 ± 0.05, and R multi = 0.66 ± 0.06 ± 0.08 
indicate a larger average suppression. The allowance of a non perfect correlation between 
the normalizations of the muons and electron data samples would certainly reduce the x 2 
value of our fit. 

The inclusion of up-going muons among the data considered, again results in evidence 
against this model. We recall the fact that the passing muons are essentially collinear 
with the parent neutrinos and that the experimental zenith angle distribution does not 
exhibit large sharp features as those predicted for example by solution (b) of [|1(J (see 
fig. 6). Therefore the relative smoothness of the passing muon data allows to exclude a 
large range of values (e, e') that correspond to few oscillations in the up-going hemisphere 
(that is 0.04 < V4e 2 + e' 2 < 0.2) and still large effective mixing e > e . 

The solution (a) of [TIJ cannot be excluded using this consideration, because its fre- 
quent oscillations do not produce sharp features given the binning of the experimental 
data, and give a constant suppression 2e 2 /(4e 2 + e' 2 ) for all zenith angle bins. The model 
has no energy dependence, and therefore this average suppression must apply to the up- 
going passing and stopping events, as well as to the up-going multi-GeV events, that 
have also a rather sharp correlation between the neutrino and muon directions. This is 
in disagreement with two features of the experimental data: (i) the passing muons have 
a suppression considerably less than both the stopping and up-going multi-GeV muons; 
(ii) the shape of the zenith angle distribution of passing muons shows evidence for a 
deformation. More quantitatively, a fit to all the data with e' = gives the parameter 
values e = 1.4 and a = 1.12 but with a total x 2 = 149 (the contribution to x 2 °f the 
throughgoing muon data being 105). 
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6 Summary and conclusions 



The survival probability Piy^ — > v^) in the case of two flavor <-> z/ T oscillations has 
a well defined dependence on the pathlength and energy of the neutrinos. In order to 
establish unambiguosly the existence of such oscillations it is necessary to study in detail 
these dependences. In the analysis of the events interacting in the detector, one can study 
a very wide range of pathlengths (10 % L < 10 4 Km) but a much smaller range of neutrino 
energies close to 1 GeV (the sub-GeV and multi-GeV samples). Therefore it is not easy 
to obtain experimental information on the dependence of the survival probability on the 
neutrino energy. In fact models where the combination L/E v (flavor oscillations), L ■ E v 
(violations of the equivalence principle) and X ~ L (flavor changing neutral currents) is 
the relevant variable for an oscillating transition probability, have been proposed as viable 
solutions of these data. Neutrino decay is also dependent on the ratio L/E v , but with a 
different functional form. 

In this study we find that flavor oscillations provide a significantly better fit to the 
sub-GeV and multi-GeV data samples than the exotic alternatives we have considered, 
however with a generous allowance for systematic uncertainties the alternative explana- 
tions can still be considered as viable. Including the upward going muons in the fit the 
alternative models are essentially ruled out. 

The upward-going muons are a set of //-induced events corresponding to much larger 
E v : for passing muons the median parent neutrino energy is approximately 100 GeV, with 
a significant contribution of neutrinos with energy as large as 1 TeV, and therefore are 
in principle a powerful handle to study the energy dependence of the neutrino survival 
probability. If flavor oscillations (where L/E u is the significant variable) are the cause of 
the suppression of sub-GeV and multi-GeV muon events, the neutrinos producing passing 
upward going muons must also oscillate, but with a smaller suppression because of their 
larger energy; moreover, for the range of Am 2 suggested by the lower energy data, one 
expects a moderate but detectable deformation of the zenith angle distribution. Both 
effects are detected. 

In the alternative exotic models we have studied here, high energy events, such as 
the passing upward-going muons, are suppressed much more (L ■ E„) than or as much 
(X ~ L) as the up-going multi-GeV events, in contrast to the experimental evidence. 

Also in the case of neutrino decay, the upward-going muon data are very poorly fitted 
by the model. Of course if neutrino have different masses (and can decay) it is natural to 
expect oscillations in combination with decay. We have not explored this possibility, but 
we can conclude that decay cannot be the dominant form of muon neutrino disappearance. 

Two results of the measurements of upward going-muons are critically important to 
allow discrimination against exotic models and in favour of usual oscillations: 

• the stopping/passing (Data/Montecarlo) double ratio for the SK upward-going 
muons |0] is r = 0.56, with a combined statistical and experimental systematic 
error of 0.07. The theoretical uncertainty in the relative normalization of the two 
sets of data has been estimated as 8% in [p5[| ; more conservatively the SK collabo- 
ration H has used 13%. Quadratically combining the more conservative estimate 
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of the theoretical uncertainty with the experimental errors, the resulting a r is 0.1. 
Therefore the suppression for the high energy passing muons, is weaker than for the 
lower energy stopping ones at more than four sigma of significance, even allowing 
for a rather large uncertainty in the theoretical prediction. This is in contrast with 
models that predict for the stopping/ passing double ratio r a value of unity (flavor 
changing neutral currents) or larger (violations of the equivalence principle). 

• The shape of the through-going upward-going muons zenith angle distribution 
shows indication of a deformation, although the no-distortion hypothesis (with free 
normalization) has a probability close to 5%. The deformation if present is a rather 
smooth one, and the distribution can be used to rule out models (such as FCNC 
with smallish e) that produce deep and marked features in the neutrino distribution 
(well mapped by the nearly collinear muons). 

The MACRO collaboration has also obtained results on upward-going muons |§ , that 
indicate the presence of an angular deformation compatible with the presence of flavor os- 
cillations (although the oscillation fit even if significantly better that the standard model 
fit is still rather poor). Preliminary results on events where upward-going muons are 
produced in (and exit from) the detector, and a second class of events that combines 
stopping upward-going muons and downward-going muons produced in the detector in- 
dicate a pattern of suppression that is only compatible with an oscillation probability that 
decreases with energy [|3~0|] . 

Also the Kamiokande collaboration p6| has measured passing upward-going muons 



with results in good agreement with Super-Kamiokande, while the Baksan collaboration 
[ fZ7f | has obtained results not in good agreement. One should also note that the 1MB col- 
laboration has in the past measured a stopping/passing ratio for upward-going muons in 
agreement with a no oscillation Montecarlo prediction |28j (see for a critical analysis). 

In conclusion, we find that the present data on atmospheric neutrinos allow to deter- 
mine some qualitative features of the functional dependence of the disappearance probabil- 
ity for muon neutrinos. This probability (smeared by resolution effects) increases with the 
pathlength L producing the up/down asymmetry that is the strongest evidence for physics 
beyond the standard model. The difference in suppression between the sub(multi)-GeV 
muon events and the higher energy through-going muons indicates that the transition 
probability decreases with energy. These results are in agreement with the predictions 
of Vp «-> v T oscillations and in contrast with several alternative exotic models. If flavor 
oscillations are indeed the mechanism for the muon neutrinos disappearance, additional 
data with more statistics and resolution (in L and E u ) should allow to study in more 
detail the oscillatory structure of the transition probability as a function of the variable 
L/E u , unambiguosly determining the physical phenomenon. It is natural to expect that 
the oscillations involve all flavors and that electron neutrinos participate in the oscillations 
(with a reduced mixing because of the Chooz limit). The resulting flavor conversions will 
have a more complex dependence on the neutrino path and energy E v \ the detection of 
these more subtle effects could become the next challenge for the experimentalists. 
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Figure 1: Ratio Data/Montecarlo for the e-like events of Super-Kamiokande. The dot-dashed lines 
are straight line fits (independent for each panel) to the data points; the dashed lines are the result of 
fitting all data (included in fig. 1, 2 and 3) with a common normalization (a = 0.84); the solid lines are 
the result of a calculation with Z/„ <— > V T oscillations with our best fit parameters for maximal mixing 
(Am 2 = 3.2 X 10—3 cV 2 and a normalization a = 1.15). 
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Figure 2: Ratio Data/Montecarlo for the /x-like events of Super-Kamiokande. See fig. 1 for a descrition 
of the lines. 
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Figure 3: Ratio Data/Montecarlo for passing and stopping upward-going muons in Super-Kamiokandc. 
See fig. 1 for a descrition of the lines. 
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Figure 4: Distributions in L/E u for the different event classes considered. In the top panel we show 
the survival probability P[ya — > for our best fit with flavor oscillations to all data (solid line) and 
the best fit with neutrino decay to the sub-GeV and multi-GeV data (dashed line). 
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Figure 5: Distributions in L ■ E u for the different event classes considered. In the top panel we show 
the survival probability P{y^ — > Vy) for our best fit to the sub-GeV and multi-GeV data. 
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Figure 6: Distributions in cos0 u for the different event classes considered. In the two top panels we 
show the survival probability P{P[i — > v^) of the two best fit points as calculated in |l6[ . 
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